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We report on measurements of the sound velocity and attenuation in a single crystal of the candidate quantum- 
spin-ice material Yb 2 Ti 2 07 as a function of temperature and magnetic field. The acoustic modes couple to 
the spins magneto-elastically and, hence, carry information about the spin correlations that sheds light on the 
intricate magnetic phase diagram of Yb 2 Ti 2 C >7 and the nature of spin dynamics in the material. Particularly, we 
find a pronounced thermal hysteresis in the acoustic data with a concomitant peak in the specific heat indicating 
a possible first-order phase transition at about 0.17 K. At low temperatures, the acoustic response to magnetic 
field saturates hinting at the development of magnetic order. Furthermore, mean-held calculations suggest that 
Yb 2 Ti 2 C >7 undergoes a first-order phase transition from a cooperative paramagnetic phase to a ferromagnet 
below T « 0.17 K. 


PACS numbers: 75.50.-y, 72.55.+S 


I. INTRODUCTION 

It has been proposed that a class of rare-earth pyrochlore 
magnets, dubbed quantum spin ice, can host Coulombic 
quantum spin-liquid 1,2 phases at low temperatures. Such 
phases support novel emergent excitations, such as photons 
and magnetic monopoles. 1-9 This has led to an active search 
for material examples exhibiting “emergent quantum electro¬ 
dynamics” in condensed-matter systems. 1,3-7,10-20 An essen¬ 
tial ingredient for the proposed rich phenomena is quantum- 
mechanical tunnelling between the macroscopically degener¬ 
ate classic spin-ice (“2-in-2-out”) ground-state spin configu¬ 
rations. 

In this context, Yb2Ti20 7 has emerged as one of the leading 
candidates for a quantum-spin-ice realization. In Yb2Ti20 7 , 
the crystal field surrounding the Yb +3 ions, that form the py¬ 
rochlore spin-network, stabilizes a well-separated S =1/2 
Kramers doublet which accounts for the low-energy mag¬ 
netic properties. 21 Recent experiments reveal that in addition 
to the dominant frustrated (classic) Ising exchange interac¬ 
tions along the local [111] directions, significant transverse 
spin exchanges are present in Yb2Ti20 7 . 10-13 The exact role 
of the transverse terms in shaping the low-energy properties 
of Yb 2 Ti20 7 is, however, far from understood. In the tem¬ 
perature window between approximately 0.24 and 1 K, recent 
experiments indicate that Yb 2 Ti 20 7 is in a correlated para¬ 
magnetic phase. 10,121 ’ However, whether this phase extends 
down to lowest temperatures giving rise to a quantum-spin- 
ice ground state is a topic of active current research. 

Notably, aside from the mentioned Coulomb spin liquid, 
the transverse terms can potentially lead to magnetic ordering 
as they lift the classical degeneracy. 10,22 Indeed, recent muon- 
spin relaxation, 23 magnetization, 24 magnetic susceptibility, 25 
and single-crystal neutron-scattering 22 measurements show 
evidence of long-range ferromagnetic spin correlations be¬ 


low T ~ 0.24 K with a reduced ordered moment. However, 
a host of other experiments have failed to detect such low- 
temperature ferromagnetic order, 26-30 although some short- 
range correlations have been observed and the paramagnetic 
state was also found to transform into a fully polarized state 
on applying a modest magnetic field. 29 Furthermore, the low- 
temperature properties of Yb 2 Ti 2 0 7 are known to be sensi¬ 
tive to small variations in stoichiometry depending on the 
oxygen content and possible Yb 3+ substitutions at Ti 4+ po¬ 
sitions (“stuffing”). 10,31,32 From the theoretical point of view 
several exotic scenarios have been suggested for the phase 
below T ~ 0.24 K including a first-order transition from a 
quantum-spin-ice state to a ferromagnet associated with the 
condensation of magnetic monopoles 10 or a low-temperature 
Coulomb ferromagnetic phase that has finite magnetization as 
well as gapless emergent photons. 4,33 

In this work, we were able to obtain a complementary in¬ 
sight into the low-temperature spin physics of Yb 2 Ti 2 0 7 by 
probing the elastic modes which couple to the spins through 
magneto-elastic coupling and, hence, can reveal information 
regarding the nature of the magnetic state, particularly the 
spin fluctuations. Indeed, previous studies showed that the 
coupling between the phonons and spins in classical spin ice 
could lead to distinct signatures in ultrasonic measurements 
which shed light on the low-temperature state and particularly 
its excitations (both in the presence and absence of an exter¬ 
nal magnetic field). 34,35 Such investigations, however, have so 
far rarely been done for quantum-spin-ice candidates. Here, 
we report on pronounced features in the acoustic properties 
of Yb 2 Ti 20 7 at low temperatures that arise due to magneto¬ 
elastic coupling. In particular, our results show evidence for 
spontaneous and field-induced phase transformations. In zero 
magnetic field, we find a softening of the sound velocity (and 
a simultaneous rise in the sound attenuation) with decreasing 
temperature due to the development of spin-spin correlations 
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in the compound. The qualitative character of the softening 
changes at around 0.17 K and shows pronounced thermal hys¬ 
teresis (Fig. 1). At low temperatures, the sound velocity (at¬ 
tenuation) shows non-monotonic behaviour as a function of 
magnetic field and a clear minimum (maximum) is observed 
at and above T ~ 0.17 K (Figs. 3 and 4). For lower tem¬ 
peratures, the behaviour is monotonic. For all cases the ve¬ 
locity and attenuation saturate for fields greater than about 
2.2 T indicating complete polarization of the magnetic mo¬ 
ments. Using a mean-field calculation, our results suggest 
that Yb 2 Ti 207 undergoes a first-order transition into a ferro¬ 
magnetic state from a cooperative paramagnet at T « 0.17 K, 
where a peak in the specific heat is also observed in confor¬ 
mity with earlier results . 10 - 22 The thermal hysteresis observed 
in the sound velocity and attenuation is a direct result of this 
first order transition. Using these insights, we draw an empir¬ 
ical phase diagram for Yt^TUCU as a function of temperature 
and magnetic field. 


II. EXPERIMENTAL DETAILS 
A. Sample preparation 

Stoichiometric quantities of Yb 203 (4N) and TiCL (4N) 
were well mixed in ethanol, dried and pre-synthesized at 
1000 °C for 24 h in air. Commercially available Yb 203 was 
pre-heated at 900 °C for 10 h in air to get rid of water. The 
pre-sintered powder was pressed into rods under 3 kbar pres¬ 
sure using a cold isostatic press (EPSI) and finally synthe¬ 
sized at 1400 °C for 24 h in oxygen flow. The crystal was 
grown at 5 mm/h and 4 bar oxygen pressure by a crucible-free 
floating-zone method with optical heating using a 4-mirror- 
type image furnace. The crystals had a typical diameter of 
5.5 - 6.0 mm and a length up to 20 mm, with yellowish 
color and well transparent in the visible light. The inves¬ 
tigated Yb 2 Ti 2 C >7 single crystals have a cubic crystal struc¬ 
ture with space group Fd3m consistent with literature 31 as 
confirmed by powder and single-crystal diffractometry. The 
microstructure of the crystal was characterized by scanning 
electron microscopy (Zeiss EVOMA15) along with compo¬ 
sitional analysis using an electron microprobe analyzer for 
semi-quantitative elemental analysis in the energy-dispersive 
x-ray (EDX) mode (X-MaxN20 detector from Oxford In¬ 
struments with AZtecEnergy Advanced acquisition and EDX 
analysis software). The composition is determined by aver¬ 
aging over several points (on the order of 10 points). The 
crystal was found homogeneous in composition and had no 
visible inclusions as confirmed by EDX analysis and optical 
microscopy. The composition of the crystal is stoichiomet¬ 
ric with respect to Yb and Ti within the typical error bars of 
an EDX experiment. The systematic error typically given for 
EDX analysis is on the order of 1-2 %. It is known that the 
long-range ferromagnetic order in Yl^TUCU is quite sensitive 
to Yb deficiency . 10 Please note that the measurement of the 
oxygen content is not reliable with EDX. 


B. Details of the measurements 


Two parallel (111) surfaces were cut and polished giv¬ 
ing a sample length of about 3.7 mm for the ultrasound- 
propagation experiments with a size of a few millimeters in 
the plane normal to the sound-propagation direction. The 
sample orientation was checked using x-ray Laue technique. 
Two acoustic modes have been studied, the longitudinal Cl, 
cl - (cn + 2 ci 2 + 4c44)/3 with k||u||[111] and the transverse 
cj mode, cj = (cn + C 44 - ci 2)/3 (k||[l 11], u±k). Here, 
k and u are the wave vector and polarization of the acoustic 
wave, respectively. The elastic constants cr are related to the 
sound velocity, vr, cr = pv\, where p is the mass density. Res¬ 
onance LiNb 03 and wide-band PVDF-film (polyvinylidene 
fluoride film) transducers glued to the sample were employed 
for sound generation and detection. The relative change of 
the sound velocity, Av/v, and the sound attenuation, A a, were 
obtained using a phase-sensitive-detection technique . 36 Cal¬ 
ibrated Ru0 2 resistor, directly attached to the sample was 
employed for thermometery. We used a dilution refrigerator 
(0.02 - 2.5 K) combined with a commercial 20 Tesla super¬ 
conducting magnet. Before each field sweep, we demagne¬ 
tized the magnet at sample temperature of ~ 2.5 K and then 
cooled the sample to the desired temperature. The measure¬ 
ments were performed under the zero-field-cooled (ZFC) con¬ 
ditions. 


Thermal expansion of Yl^TUCU was studied along the 
[ 111 ] direction by capacitive dilatometry in the temperature 
range between 60 mK and 1.5K using a tilted-plate miniatur¬ 
ized dilatometer 37 . To avoid thermal gradients the dilatome- 
ter was positioned directly into the mixing chamber of a top¬ 
loading dilution refrigerator. Vibrations were reduced by 
damping the pumping lines. All these arrangements result in 
a resolution of relative length changes of 10 9 . 


Specific heat was measured using the heat-pulse method. 
Measurements were carried out in a 3He/4He top-loading di¬ 
lution refrigerator in combination with a 16/18 T supercon¬ 
ducting magnet. The magnet was first warmed to temperatures 
above 100 K to ensure an accurate zero field measurement. 
A 332.7 pg piece of Yl^TUCU, from the same sample used 
for ultrasound measurements, was secured on a sapphire plate 
with Apiezon N-grease. The addenda, including the contribu¬ 
tion from the N-grease and the sapphire plate, was accurately 
measured and subtracted from the obtained results. A Lake 
Shore 370 AC resistance bridge was used to monitor a cali¬ 
brated thin film RuCL temperature sensor and a Yokogawa DC 
current source was used to apply a heat pulse (current ranged 
from 1-25 pA for 100-300 ms) to the 2 kOhm heater. Both the 
temperature sensor and the heater were fixed to the bottom of 
the sapphire plate which was suspended by nylon threads in 
order to provide a small heat link to the cold bath. 
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III. RESULTS AND DISCUSSION 
A. Temperature dependences 

In Fig. 1(a), we show the sound-velocity change (A v/v) as a 
function of temperature (below 0.5 K) for both the longitudi¬ 
nal and transverse polarizations. Upon cooling A v/v decreases 
untill it reaches a minimum at around 130 mK. Concomitantly, 
the sound attenuation (Arr) increases rapidly below 220 mK 
revealing a maximum at about 130 mK [Fig. 1(b)]. Below this 
temperature, both Av/v and Aa show pronounced hysterises 
down to the lowest measured temperatures. In the longitudi¬ 
nal acoustic mode the effects are much smaller than for the 
transverse mode. This may occur due to (i) indirect exchange 
interactions being dominant in Yt^TBCh, and, (ii) the smaller 
velocity of the transverse mode [which occurs in the denomi¬ 
nator of the expression for Av/v, see Eq. (15) in Ref. 38]. 

In Fig. 1(c), we plot the temperature dependence of the spe¬ 
cific heat of the single crystal where we see a clear anomaly 
at about 170 mK. This result is in agreement with earlier 
specific-heat data obtained for Yt^TUO? single crystals as 
reported by Chang et al . 10 The onset of the specific-heat 
anomaly coincides with the rapid changes in the sound atten¬ 
uation and velocity around 220 mK, whereas the maximum in 
the specific heat correlates with the inflection points for both, 
the sound velocity and attenuation. Most notably, the hystere¬ 
sis appears at temperatures lower than the specific-heat peak. 

Finally, the thermal expansion, AL/L, [Fig. 1(d)] exhibits a 
strong increase upon cooling below 200 mK without satura¬ 
tion down to lowest measured temperatures. The linear ther¬ 
mal expansion coefficient, 07 , which is giving by the temper¬ 
ature derivative of the thermal expansion, resembles closely 
the sound-velocity temperature dependence with a minimum 
at about 130 mK. Note, that AL/L is smaller than the sound- 
velocity change at these temperatures. Hence, the sound- 
velocity anomalies cannot be explained by the thermal expan¬ 
sion alone. 

In an insulating magnet, changes in the sound velocity and 
attenuation at low temperatures typically occur due to the cou¬ 
pling of the phonons to spin fluctuations. This is caused (i) 
since the sound waves modulate the position of the magnetic 
and non-magnetic ions (involved in the super-exchange) and, 
hence, modify the effective magnetic interactions, and (ii) by 
alterations of the local crystal field by changing the position 
of the ligands (nonmagnetic ions surrounding the magnetic 
one). This affects the single-ion properties of the magnetic 
ion (e.g., the effective ^-factor and the single-ion magnetic 
anisotropy). For weak sound waves and low temperatures and 
small magnetic fields (pertinent to our case), the first mecha¬ 
nism is dominant . 39 The changes in the interactions between 
magnetic ions, in turn, renormalize the sound velocity and 
attenuation. These renormalizations are then proportional to 
various spin-correlation functions 38 (see below). Hence, mea¬ 
suring the sound velocity and attenuation, both as a function of 
temperature and magnetic field, can lead to a wealth of infor¬ 
mation about the nature of the magnetic correlations present 
in the system. 

The observed features, as we shall show below, owe their 
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FIG. 1. (Color online) Temperature dependence of (a) sound veloc¬ 
ity, Av/v, and (b) sound attenuation, An, for transverse cr and lon¬ 
gitudinal cl acoustic modes, propagating along the [111] direction. 
Ultrasound frequency is 94.4 and 110.8 MHz for the longitudinal 
and transverse mode respectively, (c) Temperature dependence of 
the specific heat, (d) Temperature dependence of the thermal expan¬ 
sion, AL/L, and the linear thermal expansion coefficient, a L , along 
the [111] direction. 


origin to changes in the spin correlations as a function of tem¬ 
perature. To characterize the experimental observations, we 
generalize the calculations of Ref. 38 where it was shown that 
for cubic pyrochlores, 

Av 

—— ~ A(e sp j n DTc sp i „), ( 1 ) 

where both the magneto-elastic parameters A and D depend 
on the details of polarization and direction of the sound-wave 
propagation as well as on the sound velocity, and derivatives 
of the exchange integrals with respect to the change of the dis¬ 
tances between magnetic ions. and c sp i„ are the average 
energy per spin and the spin contribution to the specific heat, 
respectively (See Appendix A). 

Hence, as mentioned above, the acoustic signatures depend 
on the spin-correlation functions, namely two-spin and four- 
spin correlators. Approximate calculations of these correla¬ 
tors suitable for different phases then can be compared to the 
present experiments in order to obtain a better understanding 
of the magnetic phase. 
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FIG. 2. (Color online) Sound-velocity change as obtained for classi¬ 
cal spin ice within single tetrahedron approximation. The sound at¬ 
tenuation is negative of the fractional change in sound velocity (see 
text). The parameters used [with reference to Eq. (1), up to an arbi¬ 
trary shift] are A = 1, D = 1, J = 0.07. 


Further, following Ref. 39, (only keeping the q = 0 compo¬ 
nent of the magnetic susceptibility, x), the sound attenuation 
can be related to the sound-velocity change as 


6v co 2 k 2 BT X 

v v 4B 2 T 2 x 2 + w 2 ’ 


( 2 ) 


where a>k(— vk) is phonon frequency (v is the velocity of sound 
and k is the momentum) and B is a material-dependent con¬ 
stant, respectively. 

We find that the qualitative features observed in our exper¬ 
iments on Yb 2 Ti 207 , except for the hysteresis, can be repro¬ 
duced within the single-tetrahedron approximation of a classi¬ 
cal spin ice. The expressions for the specific heat and energy 
per spin are given in Appendix A. Using these expressions 
and approximating the magnetic susceptibility by its homo¬ 
geneous value, Xi for classical spin ice as that of a gas of 
monopoles 40 - 41 (see also the review article Ref. 42) giving 


a 3 k B T ’ 


(3) 


(where a is the lattice constant) we find using Eq. (2), A a oc 
-Av/v, i. e ., the sound attenuation in spin ice is proportional to 
the sound-velocity change as indeed observed experimentally. 
The temperature dependence of Av/v within our approxima¬ 
tion is plotted in Fig. 2. 

Our model can describe the qualitative behavior of the 
sound velocity, but deviates from the experiment quantita¬ 
tively. This is in accordance with the fact that we have not 
taken into account the transverse terms which are substan¬ 
tial in the actual material. Further, as noted earlier, the pres¬ 
ence of hysteresis, observed in experiments below T ~ 0.2 
K, is also not explained within this approximation. Indeed, 
hystereses have as well been observed in magnetization data 
in YbiTFOy. 24 Taking clue from these magnetization results 
along with neutron-scattering measurements, that show en¬ 
hanced ferromagnetic correlations below T ~ 0.2 K, 10,22 we 
argue that Yt^TUCF undergoes a first-order transition from 



F 0 H (T) 


FIG. 3. (Color online) Field dependence of (a) the sound velocity, 
Av/v, and (b) the sound attenuation, Aar, for the transverse acoustic 
mode ct, propagating along the [111] direction. The magnetic field 
is applied along the same direction. Arrows indicate the minimum 
positions in the sound velocity. Except for 0.02 K, data are shown 
for increasing field. The curves are arbitrarily shifted along the y 
axis for clarity. The used ultrasound frequency was 110.8 MHz. 


a cooperative paramagnet to a ferromagnet. Notably, previ¬ 
ous muon- and Mossbauer-spectroscopy measurements sug¬ 
gest that around the same temperature there is a sudden sup¬ 
pression of spin fluctuation rate from 10 3 MHz above 0.2 K 
to 1 MHz below it. 10,22,43 In the present experiments the sound 
modes have a frequency of about 100 MHz. Thus, the sound 
modes essentially see a frozen spin background and only the 
first term in Eq. (1) contributes to the fractional change in the 
sound velocity. 


B. Sound velocity and attenuation vs. magnetic field 

Figure 3 shows a field dependence of the sound velocity and 
attenuation at different temperatures for the transverse sound 
mode. At the two lowest temperatures, 20 and 120 mK, the 
sound velocity increases accompanied by a decrease in the 
sound attenuation. At 170 mK and above, a minimum in the 
sound velocity appears with a corresponding maximum in the 
attenuation. Interestingly, the changes in the acoustic proper¬ 
ties at 0.5 T are somewhat larger than the temperature-induced 
variations (Fig. 1). This can be attributed to the fact that the 
magnetic field suppresses quantum fluctuations, counterbal¬ 
ancing their effects. We also note that the hysteresis seen in 
the field sweeps is very small (only prominent in the 20 mK 
data). 

Figure 4 shows the field dependence of the longitudinal 
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FIG. 4. (Color online) Field dependence of the sound velocity, A v/v, 
for the longitudinal acoustic mode C/,, propagating along the [111] 
direction in Yb 2 Ti 2 07 . Magnetic field is applied along the same di¬ 
rection. The arrows mark position of anomalies in the low-field re¬ 
gion. The inset presents the low-temperature, low-field part of the 
H-T phase diagram extracted from our ultrasound and specific-heat 
measurements for magnetic field applied along [111] direction. Here 
FM and PM, mean ferromagnet and paramagnet. See text for details. 



FIG. 5. (Color online) (a) and (b) data and fits of the transverse¬ 
mode Av/v versus magnetic field using the asymptotic functions of 
Eqs. (4) (blue line for the paramagnetic phase) and (5) (red lines 
for the polarized phase), (c) red fit curve using a mean-field theory 
(see text) describing the experimental data at 290 mK. The fitting 
parameters are discussed in Appendix C. 


sound velocity at various temperatures. Again, similar to the 
transverse mode, a minimum (marked by arrows) is seen be¬ 
low 0.2 T for temperatures close to or above the specific-heat 
anomaly. However, the low-field anomalies become sharper 
at higher temperatures, although the corresponding anomalies 
broaden for the transverse mode (Fig. 3). 

Above the phase transition, in the paramagnetic phase, we 
can derive the field dependence of the sound velocity for small 
magnetic fields from Eq. (1) by expanding the spin correla¬ 
tions in powers of the field (see Appendix B). We obtain 

-constant + A"h 2 , (4) 

v 

where A" depends on the magnetic susceptibility and temper¬ 
ature. Similarly, for high fields (Appendix B), when the Zee- 
man term polarizes the system, we obtain 

— ~ A”'(l -D'"^)m 2 , (5) 

v \ T / 

where M ~ tanh [h/T] and x = dM/dh are the magnetization 
and the magnetic susceptibilities, respectively, h is magnetic 
field while A'" and D'" are constants. These asymptotic forms 
are applicable to general systems, beyond classical spin ice, 
and we can apply them to YbiTBOy. In Fig. 5(a), we show the 
fits to the asymptotic regions where Eqs. (4) and (5) are valid 
for the data measured at T = 290 mK shown in Fig. 3(a). The 
minimum in Av/v, corresponds to the crossover between the 
paramagnetic and polarized phase. In Fig. 5(b), we attempt to 
describe the data at T = 20 mK using Eq. (5). Clearly, at 20 
mK the response of the system is akin to that of the high-field 
behavior even for very small fields (H > 0.1 T). 

Note that, as seen from Eqs. (4) and (5), the sound- 
velocity change is typically given by the magnetization and 
magnetic susceptibility . 44 The constant value of the sound ve¬ 
locity above 2.2 T (not shown) may indicate a fully polarized 
phase where magnetization is saturated. This is consistent 
with other experiments exhibiting well developed spin waves 
about a fully polarized state in the presence of such magnetic 
fields. 1 ' However, our central observation is the fact that this 
saturation is almost complete at around 0.1 T (0.2 T) for 20 
mK (120 mK) as is evident from Fig. 3. Based on this obser¬ 
vation, we can infer that at low temperatures (T < 170 mK) 
the system is highly susceptible to a ferromagnetic phase with 
finite magnetic polarization. At higher temperature, the sys¬ 
tem enters a paramagnetic phase and requires finite fields to 
polarize. This critical field strength correlates with the posi¬ 
tions of the minima observed in Av/v marked in Fig. 3. 

We can try to interpolate between the two asymptotic de¬ 
scriptions [Eqs. (4) and (5)] within a Curie-Weiss mean-field 
theory. For that, we used Eq. (5), but with the magnetization 
being determined self-consistently, 

M - tanh [(xM + gh)/T] . (6) 

Such an interpolation gives a good fit as shown in Fig. 5(c) for 
T - 290 mK. Further, this mean-field theory predicts that the 
field required for reaching the polarized phase, characterized 
by the position of the minima (denoted by arrows in Figs. 3 
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and 4), must be proportional to temperature. Indeed, the posi- samples. A. A.Z. acknowledges the support from the Institute 
tion of the minima in A v/v grows approximately linear away for Chemistry of V. N. Karazin Kharkov National University, 
from the zero-field transition (inset of Fig. 4). 

Appendix A: Specific heat and energy of classical spin ice within 
IV. CONCLUSION single-tetrahedron approximation 


In this work, we have exploited the physics of magneto¬ 
elastic coupling to probe the magnetic behavior in the can¬ 
didate quantum-spin-ice material Yl^TUOv using ultrasonic 
elastic waves by measuring the sound velocity and attenua¬ 
tion. Our results show the development of spin correlations 
with decreasing temperatures and indicate the presence of a 
cooperative paramagnetic phase that extends down about 0.17 
K. Upon further cooling, both sound velocity and attenua¬ 
tion show non-monotonic features with a concomitant clear 
anomaly in the specific heat. The thermal expansion exhibits 
as well an anomaly at this temperature. The sound velocity 
and attenuation also show thermal hysteresis below this tem¬ 
perature. Measurements of the acoustic properties in magnetic 
fields indicate that the system is susceptible to ferromagnetic 
order for small magnetic fields above 0.17 K. Below this tem¬ 
perature, the system is best interpreted as being in a ferro¬ 
magnetic phase. Hence, we suggest that Yl^TBCU evolves 
from a cooperative paramagnet to a ferromagnet around 0.17 
K through a first order phase transition. 

However, since the sound velocity can only probe q = 0 
correlations the actual nature of the magnetic order may differ 
from that of a q = 0 ferromagnet. Further, in the light of the 
recent proposal of a possible Coulomb ferromagnetic phase 
for this compound, 4 which, along with the magnetic order, 
supports an emergent photon, we note that, unfortunately, the 
present probe is unlikely to distinguish between a regular and 
a Coulomb ferromagnet. More microscopic probes are needed 
to address this issue. It might offer an interesting avenue of 
future research. 

Finally, it may be of interest to point out a possible effect 
in ultrasonic measurements near a first order phase transition 
from a paramagnet to a ferromagnet. The ferromagnet so ob¬ 
tained will have domain walls. The elastic wave modulates 
the local exchange couplings, J,j, for the spins. 38 The sound 
waves enhance the fluctuation of the domain walls by locally 
altering the magnetic exchange energy of the bonds on which 
the domain walls lie. This then would directly affect the dy¬ 
namics in the ferromagnetic phase and the thermal hysteresis 
near the transition. Therefore, it would be of interest to see 
if such frequency-dependent effects can be observed in ultra¬ 
sound measurements. 
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For the classical spin ice, the energy per spin, s, and the 
specific heat are given by: 38 



where N s is the total number of spins. For quantum spin ice 
there are additional transverse terms. The generalization of 
above expressions is straightforward in that case. 

We can estimate the specific heat and energy for 
the classical spin-ice model within the single-tetrahedron 
approximation. 42-43 ' 46 

The partition function is given by: 


Z CC [e 2J ' T (6 + %e- 2JIT + 2e - 8y / r )f 1/2 , (A2) 


which gives 


1 <91nZ 

"arw 


— —J + 


4J (e~ 2J ' T + e~ 8J/T ) 
(3+4 e-W +e-» J IT) 


and 


C spin — 


de spin 24 J 2 e- 2J < T (l - 2 e~ 2J ' T + 3e^ T ) 


dT 


T 2 (3 + e~ 2J / T - e- AJ ! T + e~ 6J / T f 


(A3) 


(A4) 


where J is the effective exchange constant and we use units 
with kg — 1. The internal energy can be easily calculated. 


Appendix B: Derivation of the sound-velocity change in 
presence of magnetic held 

For Eq. (1), we obtain 

^ ~ A(s spin - DTc spin ) , (B1) 

where s sp j n and r v; „„ are given by Eq. (Al) and for any opera¬ 
tor O, 

(0) = ^Tr[e->™0\ 

((O 2 )) = |TV [e^O 2 ] - (O) 2 , (B2) 

with 

'H = Y j J ali sh P , - E h ' s " ( B3 > 

(if) i 
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For low fields and low temperatures (J3h < 1, fjj > 1), we 
obtain, by expanding the Zeeman term in the Hamiltonian, 


' 21,2 


P~h 


-'spin 


' s P in lh=0 ‘ 


2 N s 


V (if) 


y, h ■ s,- 


V i 


5 spm|| 1= Q + X~h~ 


(B4) 


where, M ~ tanh[/t/r] is the high-held magnetization and 
X — dM/dh is the corresponding magnetic susceptibility. This 
gives 


~ A"'(l -D'"-)m 2 , (B8) 

h V T / 

which is given in Eq. (5). 


Av 

v 


and 


'2 7,2 


P~h 


-spin ~ ^jpm | h=0 -r 

I . A 4 h 2 


Z fi ' s ' 

V (if) 


/ V i 


^P«lh=0 j2 ■ 

Adding the two contributions, we obtain, 


(B5) 


Appendix C: Fitting parameters used in Figure 5 

Here we provide the fitting parameters (up to arbitrary 
shifts) for the different panels of Fig. 5. 

a. Figure 5(a). For the blue curve at low magnetic 
fields, the fitting parameter [with reference to Eq. (4)] is 
A” as 4.2 X 10 -2 . For the red curve at the near polarized 
phase, the fitting parameters [with reference to Eq. (5)] are 
A’" as 1.0 x 10~ 4 and D'” ~ 4.5 x 1(T 2 . 


Av 

v h 


Av 

V 


+ A' 


h=0 



b. Figure 5(b). For the red curve, the fitting parameters 
( B6 ) are A”' as 8 x 1(T 5 and D'" « 3.6 x 1(T 2 . 


Thus, the magnetic-field dependence in this regime is of the 
form given by Eq. (4). 

For high fields (J3h > l,h/J> 1), the system is essentially 
in a polarized phase and in leading order of calculating the 
expectation values of spin correlations we can neglect the spin 
Hamiltonian to obtain 

&spin ~ T/ , 


c. Figure 5(c). The fitting parameters for the red curve 
[with reference to Eq. (5) where the self consistent magne¬ 
tization is obtained from Eq. (6)] are given by A’" as 0.3 x 
10~ 4 , D"' as 1.6, x « 0.1, and g as 2.5. 

Note that near the polarized phase the second term in 
Eq. (5) is subdominant for the panels (a) and (b), as expected. 
This is evident for the small values of D"' for these two fig¬ 
ures. For panel (c), D"’ is a really different constant, although 
the expression has the same form as of Eq. (5) which is why 
we have used the same notation for the coupling constants. 
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